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When polarization phenomena are considered, the usual description of the excitation state of an atomic system in terms of the 
£N) population of its energy levels is not adequate, and it is necessary to specify the population of each magnetic sublevel, as well as 
<0> die quantum interference (or coherence) that might be present between pairs of them. Whenever the magnetic sublevels of a given 
i energy level are unevenly populated and/or quantum interference between pairs of them are present, the atomic system is said to be 
,— I polarized. Atomic polarization is generally induced whenever an atomic system is excited by means of a physical process which is 
j> not spatially isotropic. 

A powerful theoretical tool that allows to describe in a very compact way the full excitation state of an atomic system is the 
/\ so-called density operator (see lFano1ll957l) . The most natural basis for defining the matrix elements of the density operator is the 
oasis of the eigenvectors of the total angular momentum \aJM), with J the total angular momentum, M its projection along the 
quantization axis, and a a set of inner quantum numbers. On this basis, the elements of the density matrix are given by 



Within the framework of the density matrix theory for the generation and transfer of polarized radiation, velocity density matrix cor- 
relations represent an important physical aspect that, however, is often neglected in practical applications by adopting the simplifying 
approximation of complete redistribution on velocity. In this paper, we present an application of the Non-LTE problem for polar- 
ized radiation taking such correlations into account through the velocity-space density matrix formalism. We consider a two-level 
atom with infinitely sharp upper and lower levels, and we derive the corresponding statistical equilibrium equations neglecting the 
contribution of velocity-changing collisions. Coupling such equations with the radiative transfer equations for polarized radiation, 
we derive a set of coupled equations for the velocity-dependent source function. This set of equations is then particularized to the 
case of a plane-parallel atmosphere. The equations presented in this paper provide a complete and solid description of the physics 
of pure Doppler redistribution, a phenomenon generally described within the framework of the redistribution matrix formalism. The 
redistribution matrix corresponding to this problem (generally referred to as Rj) is derived starting from the statistical equilibrium 
equations for the velocity-space density matrix, and from the radiative transfer equations for polarized radiation, thus showing the 
equivalence of the two approaches. 

Key words. Atomic processes - Line: formation - Polarization - Radiative Transfer - Scattering - Stars: atmospheres 



1. Introduction 



(aJM\p\a'J'M') = p(aJM,a'J'M') , (1) 

with p the density operator. The diagonal elements represent the populations of th e magnetic sublevels, the off-diagonal ele ments 
the quantum interference (or coherence) between different magnetic sublevels (see lLandi DegFInnocenti and L andolfi 2004, here- 
after LL04). For the sake of simplicity, in this paper we will only consider interference between pairs of magnetic sublevels 
pertaining to the same 7-level, which is a good approximation for the investigation of many solar spectral lines (see LL04 and 
Bell uzzi & Tru jillo Buenq |201ll for a detailed discussion on the importance of interference between different /-levels in determin- 
ing the wing polarization of multiplet lines). We will thus consider only the density matrix elements of the form p(aJM, aJM'). In 
general it is convenient to work in terms of the multipole moments of the density matrix (or spherical statistical tensors) 

p K Q {aJ)= ^(-l) y - M V2tf+l ( J M _ J M , K\p(aJM,aJM'), (2) 

MM' \ ' 



Article number, page 1 of 1 141 



which transform as irreducible tensors under a rotation of the reference system. 

Since the radiation field experienced by an atom depends, because of the Doppler effect, on its velocity, the density matrix 
will in general depend on the velocity v of the atom. Indicating with pg(aJ; v) the ^-dependent spherical statistical tensor of the 
atomic system, and with f(v) the velocity distribution function of the atoms in a given point of the plasma, a complete statistical 
description of the atom is given by the product f(v) p^(aJ; v), generally referred to as velocity-space density matrix. The need for 
the introduction of this quantity was put forward by Landi De gl' Innocentil d 1 9961) who pointed out the importance of velocity density 
matrix correlations in polarized radiative transfer. The statistical equilibrium equations for the velocity-space density matrix can be 
written in the form (see Section 13.2 of LL04) 

i [mftoj-.v)] - m(j t P>J^ o + (| [/(^^)Lc han , co , • (3) 

The first term in the right-hand side is due to processes which, as a first approximation, are not effective in changing the velocity 
of the atom: these include radiative processes (absorption and emission of photons), inelastic and superelastic collisions with 
electrons, and depolarizing collisions with neutral hydrogen atomsQ The second term, which can be regarded as a generalization 
of the Boltzmann term which is met in the kinetic theory of gases, is due to collisions that are able to modify the velocity of the 
atom (velocity-changing collisions). Such collisions, which are characterized by very small impact parameters and by rather large 
exchange of kinetic energy, generally induce transitions between different energy levels, and thus affect the atomic density matrix. 
Because of the generalized Boltzmann term, Equation Q is extremely complicated, and two different approximations are generally 
introduced. 

The first one consists in neglecting the generalized Boltzmann term (velocity-coherence approximation). As discussed in Section 
13.2 of LL04, this is a good approximation in the outer layers of a stellar atmosphere, where the number density of perturb ers 
responsible for velocity-changing collisions (typically hydrogen atoms or ions) is sufficiently low, and this kind of collisions are 
indeed negligible. Under this approximation, the ^-dependence of the density matrix is only due to the fact that atoms moving with 
different velocities may experience, because of the Doppler effect, different radiation fields. 

The second approximation is to assume that velocity-changing collisions are so efficient in reshuffling the atomic velocities that 
any velocity density matrix correlation is lost (complete redistribution on velocity approximation). In this case, the velocity-space 
density matrix is given by f(v) pg(aJ), the density matrix being independent of v. Since velocity-changing collisions also contribute 
to depolarize the atomic system, when the complete redistribution on velocity approximation is justified, polarization phenomena 
are generally negligible. For this reason, it is customary to consider an intermediate approach which consists in neglecting the 
generalized Boltzmann term, still assuming a velocity-independent density matrix. 

The general problem of interpreting the spectropolarimetric profiles of lines formed in an optically thick plasma, such as a stellar 
atmosphere, requires the self-consistent solution of the statistical equilibrium equations and of t he radiative transfer eq uations, taking 
into account polarization phenomena, both in the atomic system and in the radiation field (see Truiillo Bueno 2003). This problem 
has been referred to as the Non-LTE problem of the 2 nd kind (see LL04), so to distinguish it from the usual Non-LTE problem 
where polarization phenomena are neglected. A detailed discussion of the general Non-LTE problem of the 2 nd kind, under the 
approximation of complete redistribution on velocity, is presented in Chapter 14 of LL04. 

In this paper, we present an application of the Non-LTE problem of the 2 nd kind under the velocity-coherence approximation 
previously discussed. We consider the basic case of a two-level atom with infinitely sharp upper and lower levels. Starting from 
the statistical equilibrium equations for the velocity-space density matrix (Sect. 3), and from the radiative transfer equations for 
polarized radiation (Sect. 4), we derive a set of coupled equations for the velocity-dependent source function (Sect. 5). This set of 
coupled equations is then specified to the particular case of a Maxwellian distribution of velocities, and it is finally applied to the 
particular case of a plane-parallel atmosphere (Sect. 6). From the same equations we also derive the redistribution phase-matrix 
corresponding to this physical problem (generally referred to as Rj), thus showing the equivalence of this latter, widely applied 
approach to the one described in this paper (Sect. 7). 

The more realistic case of partial frequency redistribution, assuming a two-level model atom with infinitely-sharp and unpo- 
larized lower level, and with a naturally and/or collisionally broad en ed upper level, ha s been considered by several authors (e.g., 
lOmont et al.lll972Lll973l:lDomke & Huben\ll98^lBommierill997aK[Sampoornal20T2l) . providing expressions for the R n and R m 
redistribution matrices. In this paper we show, for the academic case of a two-level atom with infinitely-sharp upper and lower 
levels, how the phenomenon of pure Doppler redistribution can be rigorously described through the velocity-space density matrix 
formalism. As we will see, this formalism allows to describe the physics of the atom-photon interaction in a very transparent way 
(e.g., it allows to clearly identify correlations between atoms located at different points of the plasma), and it highlights the un- 
derlying approximations (e.g., the neglect of velocity-changing collisions). The equations that are obtained are very general, and 
can be applied to arbitrary velocity distributions. Moreover, they are able to account for the presence of atomic polarization in the 
lower level of the considered transition, and they are perfectly suitable for a generalization to the multilevel case. The application 
of the velocity-space density matrix formalism for the description of a two-level atom with broadened upper level is presently under 
investigation, and will not be discussed in this paper. 



1 Depolarizing collisions are elastic collisions due to long-range interactions, and are thus ineffective in changing appreciably the velocity of the 
atom. 
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2. Formulation of the problem: hypotheses and approximations 

We consider a two-level atom without hyperfine structure and with unpolarized lower level. Consistently with the last assumption, 
we suppose that the radiation field incident on the atom is weak, in the sense that the average number of photons per mode, «, is 
much smaller than unity, which justifies to neglect stimulated emission. 

We suppose that a collection of such atoms is distributed within a static medium of arbitrary shape. In this medium the atoms 
interact with a magnetic field, B, and with a population of colliding particles having a Maxwellian distribution of velocities charac- 
terized by the temperature T. No restriction is made on the spatial variation within the medium of the temperature T of the colliders, 
of the densities of the atoms and colliders, and of the magnetic field vector B. We suppose that the magnetic field is weak (in the 
sense that the associated Larmor frequency vl is much smaller than the frequency width Avp of the absorption profile) and that the 
inverse lifetime of the upper level, y„, is also much smaller than Avp, so that the flat-spectrum approximation is satisfied^ 

The atoms are characterized by an arbitrary velocity distribution f(v) that, for the sake of simplicity, we assume to be constant 
throughout the medium. Taking velocity density matrix correlations into account, at any point P of the medium, of coordinate 
x, the atom is thus described by the velocity-space density matrix f(v) [pg(aJ; i>)] , where (aJ) = (a u J u ) for the upper level and 
(aJ) = (aeJe) for the lower level. Neglecting velocity-changing collisions (i.e. assuming the velocity-coherence approximation 
described in Sect. 1), the velocity-space density matrix evolves with time according to the equation 

i {m \ftoj-, 4) = m [j t W<*j-, ^ , (4) 

which is solved by 

4/(»)=0, (5) 
df 

^[pli^v)} x =[^[p K Q {aJ-,v)\] . (6) 



o 



As previously pointed out, the term in the right-hand-side of Eq. (0 contains the "ordinary" processes due to the atom-radiation 
interaction and to collisions. The value of [pg(ar./; u)J can thus be found by solving the statistical equilibrium equations presented 
in LL04, taking properly into account the explicit dependence of the radiative and collisional rates on the velocity of the atom. The 
velocity distribution f(v), on the other hand, remains undetermined, and can only be established by means of different physical 
considerations. In many cases, it can simply be assumed to be a Maxwellian, possibly centered at a non-zero velocity, like in the 
case of the solar wind. In the next Sections, if not explicitly specified, we consider an arbitrary velocity distribution f(v). 

Finally, we neglect the broadening effect of elastic collisions, so that both the upper and lower level of the atom can be considered 
infinitely sharp. Consequently, we also neglect any frequency redistribution effect due to elastic collisions. In terms of scattering 
processes, our model thus allows us to describe coherent scattering in the atom rest frame, with purely Doppler redistribut ion in the 
observ er frame. Within the framework of the redistribution matrix formalism, following the terminology introduced by iHummerl 
(fT%l . this kind of process, in the case of a Maxwellian velocity distribution, is described by the R[ redistribution function. 



3. The Statistical Equilibrium Equations 

Referring to the geometry of figure 1, in the 'fixed' (or laboratory) reference system S, the statistical equilibrium equation for the 
multipole moments of the upper level is the same as Eq. (14.2) of LL04, with the only difference that the radiative and collisional 
rates now depend explicitly on the velocity v of the atom. We thus have 

— [p£(ar B /«;»)] J =-2itivLg a ,j u Yj'Kqq' [pQ,(a u Ju\v)\ x 
at 'o Q 

+ ^ T A (a u J u KQ,a e J ( K'Q') [p%,(a e J e ;v)] % 

K'Q' 

- [ma u J u KQK'Q') + R s (a u J u KQK'Q')] [p K Q ,(a u J u - v)\ 

K'Q' 

+ J^j- + + | Cf\a u J u , a c J t ) [p K Q (ac h \ «)] x 

- [cf(a ( Je, a u J u ) + D {K) (a u J u )\ [p K Q {aJ u - v)] % , (7) 

where all the rates are evaluated at point x, vl is the Larmor frequency at the same point, and g a „j„ is the Lande factor of the upper 
level. The kernel "Kqq, is given by (see Eq. (7.79) of LL04) 

Kq ^ Z Z&eC**)* Q" ®q»q'(Rb) , (8) 

Q" 

2 The applicability of the flat-spectrum approximation follows from the two inequalities y u <k Av P and v L <k Av P . The latter obviously implies 
an upper limit on the magnetic field intensities that can be handled by this formalism. 



Article number, page 3 of 1 141 



Fig. 1. At each point P of the medium, the magnetic field vector B is specified by the angles 6 B mi&xb, defined in the fixed reference system 
E = (xyz). The polarization unit vector e a (£2) specifies the reference direction for positive Q of the radiation flowing through P in the direction SI, 
specified by the angles 6 andx in the fixed reference system. 



where Dqq,(R) are rotation matrices, and Rg is the rotation that carries the local 'magnetic' reference system (having the z-axis 
aligned with the magnetic field) into the 'fixed' reference system E. In terms of Euler angles one simply has (see Fig. 1 for the 
definition of the angles) 

Rb = (-7b, -8b, ~Xb) , 

where jb is an arbitrary angle that can be set to zero. The main properties and the explicit expressions of the components of 'Kqq, 
are given in App. 19 of LL04. The assumptions that we have introduced yield two basic simplifications in Eq. ||7}: 

- because stimulation effects are neglected, the relaxation rate Rs is zero; 

- because lower-level polarization is neglected, the statistical tensors of the lower level reduce to pg(a(J(; v) = p^iaeJf, v) 5kq Sqq. 

This implies that the only radiative rate needed in Eq. Q is the absorption rate, Tp,(a u J u KQ, a(J( 0). Taking into account that the 
upper and lower levels of the atom are infinitely sharp, such rate, in the fixed reference frame, is given by (cf. Eq. (10.9) of LL04) 

T A (ar»/«*fi,a^00)= ^{2J e + 1) B(a t J t -» a u J u )(-\) l+] ' +] "+ Q j ^ ^ J J [/ fi (vo)] cf , (9) 

where |/5g(vo)] f is the radiation field tensor calculated in the atom rest frame, or comoving frame (it describes the radiation field 
as "seen" by the atom, and it actually depends on its velocity v in the fixed reference frame), and where vo is the frequency of 
the transition between the upper and lower level of the atom. In the comoving frame, the radiation field tensor is given by (see 
Eq. (5.157) of LL04) 

3 

[j- Q (va)] c f . = 9 Z r5 Q (i, n) [Si(v , n)] c . f . , (io) 

J 71 i=0 

where [5,(vo, ^)] c .f. are the Stokes parameters of the radiation propagating along the direction Q at the frequency vq, as defined in 
the comoving frame. In the limit v/c <k 1 aberration and further relativistic effects can be neglected and the radiation field tensor in 
the comoving frame can be calculated from the Stokes parameters in the fixed reference frame just taking into account the Doppler 
effect evaluated to first order in v/c. This brings to the following expression 



!=0 



where now 5, are the Stokes parameters measured in the fixed frame. 

As far as the the inelastic and superelastic collisional rates, CJ (a u J u ,a(Ji), and Cf\a(J(,a u J u ) are concerned, we suppose 
that they do not depend on the velocity v of the atom undergoing the collision. This is a very good approximation, since such 
collisions are due to electrons, which travel with velocities much larger (typically by two orders of magnitude) than the velocity of 
the atom. Concerning depolarizing collisions (elastic collisions), basically due to neutral perturbers such as hydrogen atoms, the 
approximation can be more questionable. For this reason we will, from now on, denote such rates with the symbol D (K \a u J u ; v). 
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Taking into account these remarks, and recalling the expression of the radiative rate Re (cf. Eq. (7.14e) of LL04), which 
obviously does not depend on the velocity of the atom, being due to spontaneous de-excitation processes, we can rewrite Eq. (0 in 
the form 

- [A(a u J u -> aeJe) + cf{a e J e ,a u J u ) + D (K) (a u J u ;v)] [p*(a u J u ;v)] x 

+ Jf^y [ B{atJl ^ a " Ju) w f,i (_1)G feMcf. + 6 ko $Q0 Cf\a u J u , atJe)] [p° (a e J e , v)]^ , (12) 

where the symbol Wj 3 is given by Eq. (10.1 1) of LL04. 

Since the colliding particles have a Maxwellian distribution of velocities, we can apply the Einstein-Milne relation to connect 
the collisional rates due to inelastic and superelastic collisions (cf. Eq. (10.49) of LL04). Next we divide both members by 
A(a u J u — > a/J() and introduce the usual notations (cf. Eqs. (10.51) and (10.28) of LL04) 

, = Cf{a e J ( ,a u J u ) 5 ( K)(v)= D {K \a u Ju\v) g = 2™l g ai ,j„ 
A(a u J u -> a c J c ) ' " A(a u J„ -> a c J c ) ' " A(a u J u -> a t J c ) ' 

Recalling the relations between the Einstein coefficients (Eqs. (7.8) of LL04), we obtain, for stationary situations 

[l + e + 5f\v )] [p K Q {a u J u - v)] t + iff, J] K K QQ , [p K Q , (a u J u ; v)] s = 

Q 

V5tT [ w « (_1)G [/ -e (vo)]cf - + 6ko 6 e° eB ^ v A [p>( J r> »)], > (14) 



2hv\ V 27, 



where 



B7(vo)= _o exp (__^) (15 ) 

is the Planck function in the Wien limit (consistently with the fact that stimulated emission is neglected). 

In view of the following applications, it is convenient to rewrite the previous equation for the time evolution of the density matrix 
in a more compact form, by introducing suitable 'source functions' for the different statistical tensors. Defining 

K 2hvl FTJTTx [pfrcuMl ntr . 

o Q (v,x) = —r- - — - -j— T- , (16) 

*D° (a e J c ;v)] x 



2J U + 1 
such equation becomes 

[l + e + 8f \v )\ S K Q (v, x ) + i H u J] <K* Q , S K Q ,{v, x ) = wg> (-1) G [/5 fi (v )] cf + 6 K0 eB T (v ) . (17) 



Q 



The quantities Sq(v, x) are the obvious generalization of the irreducible components of the two-level atom source function, intro- 
duced in LL04. Now we also have an explicit dependence on v. This is because atoms having different velocities may have, in 
general, different source functions since, even at the same point in the medium, they experience, due to the Doppler effect, different 
radiation fields. 



4. The radiative transfer equation 

We consider now the radiative transfer equation. From Eq. (6.83) of LL04 we have, neglecting stimulated emission 
, 3 

— S,(v, fl) = - Yj k u S J (v > a) + ei (i = 0, . . . , 3) , (18) 

7=0 

where Sj(v, £1) are the Stokes parameters of the radiation flowing through point x in the direction £1, defined with respect to the unit 
vectors e fl (fi), e^(Q) of Fig. 1, Kf- is the absorption matrix, and e, are the emission coefficients in the four Stokes parameters. The 
explicit expressions of the radiative transfer coefficients for the case we are concerned with can be derived by a simple generalization 
of the results contained in Section 14.2 of LL04. Since we have assumed that the lower level is unpolarized, and that the Zeeman 
splitting is negligible with respect to Avd, the absorption matrix Kfj is proportional to the identity matrix, i.e., it is of the form 

^ = 77 A (v.«H7- (19) 
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Fig. 2. A ray with direction SI enters the medium at point x , where its Stokes parameters are S. (v, SI). At point x, the Stokes parameters are 
given by Eq. ( 124b . 



Due to the Doppler effect, an atom with velocity v absorbs radiation propagating in direction £1 only at the frequency vo(l + v ■ fi/c) 
(we assumed that the upper and lower levels are infinitely sharp). The absorption coefficient t]q(v, £1) is thus given by 

77 A (v,") = fc L A (x)p(v,£i), (20) 
with 

p(v, £1) = J d 3 vf(v) 6 jv + v V —^- - vj , (21) 
where 6 is the Dirac-delta. The quantity £ A (x) is the frequency-integrated absorption coefficient of the line, given by 

£ A (x) = ^ Nf(x) B{a e Jc -> aj u ) , (22) 

with Nt(x) the number density of atoms in the lower level at point x. 

The expression for the emission coefficient in the four Stokes parameters is obtained from Eq. (7.16e) of LL04 (where the 
velocity-independent density matrix has now to be substituted with the velocity-space density matrix, and where, consistently with 
our assumptions, the profile </>(vo - v) is now a Dirac-delta). Considering that, due to the Doppler effect, the atoms having velocity 
v emit, along the direction £1, at the frequency vo(l + o ■ fl/c), we have 

Si {v, £1) = £ A (x) f d 3 v f(v) $ (vo + v — - v] J] ^ ( T*(i, £1) S K Q {v, x) . (23) 

J V ° I KQ 

Obviously, the preceding expressions imply that the only contribution to the opacity and emissivity of the medium comes from 
transitions between the two levels of the model atom. The case where a source of continuum opacity (and emissivity) is also present 
is formally more complicated and will not be treated here. 

The radiative transfer equation can be formally solved. Referring to Fig. 2, and using Eq. (8.18) of LL04, the Stokes parameters 
at point x of the radiation at frequency v, flowing along the direction £1 can be expressed in the form 

Si(v, £1) = f ds' f d 3 v f(v) 5 (v„ + vo — - v) £ A (x ') T w ® T K Q (i, £1) S K Q (v, x ') + 5, (b) (v, O) , (24) 

where 5 ; (b) (v, £1) is the Stokes vector of the radiation entering the medium at point Xq along the direction £1, s' is the coordinate of 
x ' reckoned along £1 (.?' = |x ' - xo|), and r v (x, x ') is the optical depth at frequency v between points x and x ', 

t v (x,x') = j ds"j] ( f(v,£l) = j di"/t A (x") J d 3 vf(v)6\v + v ^ - vj . (25) 

It is now possible to find the expression for the radiation field tensor at point x. Substituting Eq. (f24l > into Eq. (fTTb . we obtain two 
contributions and we can thus write 

[4(vo)] cf = [4(v )] I + [Jq(vo)] e , (26) 
where the 'internal' part [■^q(vo)] i is given by 

K(4= <f ^ J] fds' f dV/C'Wvo— -n) ^ A (x') e-"^) V wZ]T%{i,£l)S%{v',x'), (27) 

J 47r i=0 Jx J \ c j K , Q , 



where 



v ■ £1 

V = vo + v , (28) 

c 
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and the 'external' part [-^(vo)]^ originating from the boundary conditions, by 

[4(vo)] E = j> ^ J) T%(i, «) e -^<**>> S,. ( V, O) . (29) 

Equation (l27l i can be cast in a simpler form by changing the double integral in dQ and ds' into a volume integral. Since 

dV = (x-x') 2 dQds' , (30) 

we get 



5. Coupled equations for the velocity dependent irreducible components of the source function 

We can now substitute the expression of the radiation field tensor at point x into the statistical equilibrium equation. From Eqs. (flTt , 
( 1261 ). and ( TJTT ) we obtain 



[l + e + tf <*>(»)] X) + iH u J] ^QG' 5 2'( ' x ) = &° <%> £B t(vo) + wfX (-D e [j- Q (?oj\ B + 

a 

+ J d ' x ' 4 ^l } 2 fd 3 v'f(v')J]G mK , Q ,(v,x-,v',x')S%,(v',x'), (32) 



where 



G KQ jc>q>(v,x;v',x') = e~ 



6\vo- 



(33) 



i=0 



The quantities Gkq,k' Q'(v,x;v',x') appearing in this equation represent a factor (having the dimension of the inverse of a frequency) 
which weights the amount of coupling between the statistical tensor pg of the atoms having velocity v at point x and the statistical 

tensor pg, of the atoms having velocity v' at point x'. They are a generalization of similar quantities introduced in LL04 and can 
be referred to as velocity dependent multipole coupling coefficients. The Dirac-delta appearing in their expression is responsible for 
the fact that two such multipoles can be coupled only if the velocity difference, v - u' is perpendicular to the unit vector Q which 
specifies the direction x-x'. 

Equation d32l is a system of linear, non-homogeneous, integral equations in the unknowns Sq{v,x), the velocity dependent 
irreducible components of the source function, which can in principle be solved once the properties of the medium and the boundary 
conditions are specified. When the values of these components are known at each point, the Stokes parameters of the radiation 
emerging from the medium can be computed by applying Eq. (l24l i. It should be remarked that, owing to a property of the tensor 
Tq , whose proof can be found in App. 20 of LL04, it can be shown that Eq. ( 1321 ) decouples in two different sets of equations 

involving, respectively, the components with K = 0, 2 and those with K = 1. In the latter set, the only source term is [/l e (vo)]E, 
which vanishes unless the boundary radiation field has some contribution arising from circular polarization. Excluding this case of 
limited interest, all the components S 1 q(v, x ) are everywhere zero in the medium. 

For further developments it is however more practical to rewrite Eq. ( l32t in an alternative form by substitution of Eq. (|33T > and 
by introducing the compact symbol Tkq,k'Q'(^)> whose main properties are collected in App. 20 of LL04. Its definition, that we 
recall here for completeness, is the following 

3 

T KQ ,K- a {m = Yj (-D G T -Q& ") T Q'(>> ") • ( 34 ) 

i=0 

With these transformations we obtain 

[l + € + 6f(v)]s^(v,x) + iH u Y i K%Q,S%{v,x) = 6 K0 6 Q peB T {v () ) + ^(-l) 2 [/ 2 (v )] B + 

Q' 

+ f dV 4n(x- If I dV/(»v(vo^-fi).Z ■ 05) 

This is a very general set of coupled equations for the velocity-dependent source function that in principle can be solved numerically 
by suitable discretization of Sg(v, x) over the velocity-space and the physical-space. 
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We now assume that the velocity distribution /(») is a Maxwellian characterized by the thermal velocity v t (possibly containing 
the contribution of microturbulent velocities) 

m^^—e-" 1 !"-. 06) 

V, 7T ' 

For fixed x and x' (which implies a fixed direction Q joining x' with x), and for a fixed velocity v, we introduce a right-handed, 
Cartesian coordinate system (e a , eb, e c ) in the velocity space. The unit vector e c is directed along £2 (coinciding with it), while the 
other two vectors are perpendicular to Q and for the rest arbitrary. In this system we obviously have 

v = u. d e a + v b eb + v c e c , v ' = v' a e a + v' h eb + v' c e c . (37) 
With these position, we can perform the following formal substitution in the integral in d 3 v' 

'di>: f do' rdv' e -W 2+l, b 2 +< 2 )/"r 



r . , ( v'-v \ rdv' rdv' rdv' H +<- m i v '-v c \ 



(38) 



Due to the presence of the Dirac's delta, the integral in dv' c is immediately performed, and the last term in the right-hand side of 
Eq. ( f35l > acquires the form 

where the vector «„' in the argument of has components (v' a , v' b , v c ), and where the Doppler width Avd is given by 

Av D = vo- . (40) 
c 

We also observe that under the assumption that the velocity distribution is a Maxwellian characterized by the thermal velocity v t , 
the absorption coefficient (see Eq. d20i i) does not depend any longer on the propagation direction of the radiation (fi), and it is given 
by 

77 A (v) = fc L A (x)p M (v-vo), (41) 
with 

Pm(v - vo) = - 7 J— e-^ 2 ^ . (42) 
Under the same assumption, the optical depth t v (x,x') (see Eq. d25l )) is given by 

r v (x, x ') = Pm(v - vo) f ds"k£(x ") . (43) 

6. Equations for a plane-parallel semi-infinite stellar atmosphere 

We now consider the particular case of a plane-parallel, semi-infinite stellar atmosphere. In this case all the physical quantities of 
the medium depend on a single coordinate, the height in the atmosphere, that we assume as the z-axis of our fixed reference system 
of Fig. 1. As a consequence, the irreducible components of the velocity-dependent source function only vary with the height z. We 
keep assuming that the velocity distribution is Maxwellian and we introduce the line optical depth through the equation 

kHz) 

d/ L = -^— dz. (44) 
Av D 

Assuming that the stellar atmosphere is not illuminated by external sources of radiation, Eq. (l35t takes the form 

[l + e + 5f\v )] S K Q {v, f L ) + i H u J] <Kl Q S K Q ,(v, f L ) = d K0 5 QQ eB T (v ) + 

Q 

pea poo poo i p j / p -(j 1 ' 2 + v ' 2 +v ^)/ l f 

where the indices K and K' are restricted to the values and 2. 

The integral over x' and y' can be performed by introducing cylindrical coordinates and following a procedure similar to the one 
developed in App. 21 of LL04. Referring to the geometry of Fig. 3 (which represents the case ti > or z' < z), and introducing the 
cylindrical coordinates r and a of the point P', one has 

dx' dy' — r dr da . (46) 
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Fig. 3. In the reference system £ we consider a fixed point P located at height z corresponding to line optical depth t L . The point P' lies on the 
plane T parallel to the x-y plane; its height is z' corresponding to line optical depth t' L . The cylindrical coordinates (r, a) define the position of P' 
in the plane T. The angles 6 and^ specify the direction £1. 

On the other hand 

(x-x'f = (z-z') 2 + r 2 , (47) 
and 

r=(z-z')tm6, Ar=^—^A8. (48) 

cos l 8 

From these relations we get 
Ax' Ay' 



(x-x') 2 



= tan 8 A8Aa , (49) 



and since \ — a + n, the double integral over x' and y' can be transformed into an integral over the angles 8 and \ specifying the 
direction Q. In Eq. d45b one can then perform the formal substitution 

Ax' \ Ay'— —r > — I d* d0tan0. (50) 

oo J -co 4tt(x-x') 2 4ttJ Jo 

Moreover, introducing the reduced frequency distance from line center through the usual expression 

(51) 

Av D 



and the normalized profile 
1 



^) = ^e-f\ (52) 



one has 

Ty(X,X ) = , (53) 

cos 8 

with (recalling Eq. (l28l l) 
Avd Avd c 

Taking into account these transformations, and performing a similar analysis for the case t' L < fL, Eq. ( l45l l becomes 
[l + e + 6f\v )]S K Q (v, tO + i H u J] 'Kqq S k & (v, f L ) = S KQ 6 Q0 e B T (v )+ 

a 

X°° 1 r^ n r Hi/ r An' <»-K 2 + i >h 2+i ' 2 )/ 1 ' 2 « . 

d^lf *xj d.ltan.1 e -^O^s 1^ ^ ^ w % W ff rWQ)gg(..',£) » W 

0\ t t K'Q' 

the interval (8 l , 8 2 ) being (0, n[2) if t' L > f L and (tt/2, n) if f L < f L . 
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It is now necessary to specify the unit vectors e a and eb that have been left undefined. For a given direction £2 we set 

e a = eg = cos 8 cos^i + cos 8 sin^j - sin 8k , 
e b = e* 



sin^i + cos^j 



(56) 



e c = e r = SI = sin6> cos^i + sin# sin^-j + cos 0k . 

Given the velocity components v' a and v' b , the velocity is thus given by 

K - v 's. e a + u b e b + y c e c = [(cos # u a + sin 8 v c ) cosy - sm Jf v'h\ • 
+ [(cos # u a + sin y c ) sin^f + cos _y u b J j 
- [sin 8 v' a - cos 8 v c ] k . 

This equation allows us to find the modulus, v' t , polar angle, 8*, and azimuth,^,, of the velocity v' t in terms of u' a , v' b , and v c . One 
gets 

v' t cos 8* — - sin v' a + cos 8 u c , 

sin 6, cosy* = (cos i£ + sin u c ) cos^- - sin^- u' b , (57) 
v' t sin 0* sin^» = (cos 8 v' a + sin 8 v c ) siny + co&x v' b , 

the combination of the last two equations giving 

v' t sin ft, e 1 ** = (cosSv^ + sin0u c + iv' b )e' x , (58) 



cos 8 v' + sin 6 v c - l v. 

,-ur. - ? b 



J(cos 8v' a + sin u c ) 2 + i> b 2 

The previous equations can also be inverted to give 

u' a - u' t sin 8 t cos 8 cos(x* — X) — v'„ cos 8* sin 8 , 
u b = v t sinfl. sin(x* . 

- y * sm #* sin # cos(^„ — x) + v' t cos ft, cos 8 . 



(59) 



(60) 



We now consider the simplified case of cylindrical symmetry. This implies the absence of a deterministic magnetic field (the case 
of a turbulent magnetic feld can be ha ndled with a slight modification of the formalism; see Sect. 14.2 of LL04, and Appendix A of 
iTruiillo Bueno & Manso Sainzll 19991) . In a cylindrically symmetric environment the velocity-dependent irreducible components of 
the source function, Sq(v, tO, have a dependence on the azimuth^ of the velocity of the form 



S K Q (v, tO = S%v, G v ,Xv, tO = S K Q (v, 8 D , 0, tO e" 



iQx,, 



(61) 



where Sg(v, 0, tO is the value of the source function corresponding to Xv = 0. Taking into account this property, we can limit 
ourselves to find the coupled equations for the quantities Sq{v, 8„, 0, tO, thus finding 



[l +e + 8f ) {v)]sl(v,6 v ,Q,tO = SkoSqo cB t (v ) 

+ i^rj *rJL«i»«i 



,-«-«L)*<f)/cose 



r dv^ rdu' b e -K 2 H 2+ >i)/»i 



(62) 



K'Q' 



where, according to Eq. 



cos 8 v, + sin 8 v r - l or 



tJ(cos 8v'. d + sin u c ) 2 + i> b 2 

We remind that the quantity v c appearing in this equation is defined by 
v c = v ■ SI = v (cos 8 V cos 8 + sin 8 V sin flcos^-) . 
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From Eq. (l63l l. one can notice that the integrand in Eq. (l62t is an odd function of vL, When integrating over v! from — oo to oo, one 
can thus perform the following substitutions 



cos 8 v' a + sin 9 v c - i v' b 



yj(cos 9v' R + sin 8 v c ) 2 + v. 



cos 8 v' + sin 9 u c - l v. 



cos 8 v' + sin 8 v c 



^(cos 8 v' R + sin 8 v c ) 2 + v' 2 



±2 



yj(cos 9 v' a + sin 8 v c ) 2 + u' b 2 

(cos 9 v' a + sin 9 v c ) 2 - v' b 
(cos 9 + sin 8 v c ) 2 + v' b 



(65) 



Finally, one can notice that the structure of Eq. (162j i is such to be consistent with our hypothesis on the behavior of the source 
function with the azimuth of the velocity. This is easily proven by considering the fact that the quantity Ykq.k'Q'(Q) depends on the 
angles through an exponential of the form exp[i(2' - Q) X ]- 

It is convenient to perform a change of variables in the integral appearing in Eq. (l62l >. Once the values of v and 9„ (the velocity 
vector appearing as the argument of the density matrix element for which we write the statistical equilibrium equation) and the 
direction Q (through the angles 9 and^f) are specified, the velocity component v c is fixed. We can then transform the double integral 
in the variables (dv' a , dv' b ) in a double integral over the variables 9*). For this we have to consider the formal transformation 



J\ 

where J' is the determinant of the Jacobian of the transformation, namely 
dv' dv' ^ 



(66) 



J" = det 



d(cosft«) d(cos&,) 



dvL 



dvl 



(67) 



On the other hand, from the equations relating v' t and 9 t with v' a and v' b , we have 



dv'a v' t ' 



dv'„ 
dv' h 



d(cosft») -sin# 



dv' 



d(cos6» t ) 
dv' h 



= 0, 



(68) 



so that we obtain 

,,'2 



1 



v: 



sin 9 v. 



(69) 



This equation, together with the relationships previously developed among the components of the velocity, allow to rewrite Eq. 
in the form 



J v t J v t J v t v t J 



1 



sin 6» | sinOf * ~x)\ ' 



where 



\sin(x* -X)\ 



1 



- u« cos 9„ cos I 
v'„ sin 8, sin 9 



(70) 



(71) 



Substituting this result into Eq. (1621 1. and inverting the order of the integrals, it is possible to rewrite the same equation in a different 
form. Performing the following change of notations on the integration variables: v' t — > v', 9* — > 9' v , one gets 



dt' L — - d9 , v YM%Xv,9 D J L ,v',9' v j' L )S%,(v\9' v ,0,t' L 

J ft ft J £7q, 

where the kernel, ZR^,(y, 8 V , t L , v', 9' v , t' L ) is given by 



e v , t h , v', t'0 = l^- 2 £dx £*e ^ e^*™- ^ w % r W «) ■ 



e -iQ'x, 



sin(x« - X )\ 



(72) 



(73) 



and can be more conveniently expressed in the form 

-In 

cos ( 



i i r 2n r 02 i 

M™(v, Bo, t L , v', 8' u , t' L ) = — — d x \ d8- 

uu 4n n 3 ' 2 J J g | cos 



e v 'l 



(^-fLj^O/cose p -(v'M\ A ,(K) (?) 



|sin(Y, -X)\ 
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(74) 

The quantities cos(x* - x) an d sinCf* _ x) appearing (implicitly or explicitly) in this equation are related to the different variables 
through the equations 



cos(y* -x) = 



v c - v cos ft;, cos ( 
u' sin ft' sin ft 



sin(y. -X) = ± A 1 - 



v c - v' COS ft' COS ( 



v' sin ft' sin ft 



(75) 



(76) 



Concerning the + sign appearing in this last expression, it is important to note that it is connected with the sign of the velocity 
component, v' b , that has now disappeared from the equations, due to change of variables that we have performed. Equation d74l has 
then, more properly, to be written as 



cos ft 

KQX'Q 



e -(£-, L W)/cos0 e -V/v t ? w (K) W (T) X 



1 



e -i<I(Xi'-X) +e ~iQ'(X. -X) 



(77) 



where x* } and^i +) are the values of^* corresponding, respectively, to negative or positive vi. With easy algebra we find 



e -iQ'(x',- ] -x) + Q-iQ'ix^-x) 



e -iQ'(xV-x) + t -iQ'(x iM -x) 



e -iQ'<xt ) -X) + e -iQ'(x ( *'-X) 



= 2 , for Q' = , 

= 2cos0t,-^) , fore'=±l , 

= 2[2cos 2 (^- ;r )-l] , fore' = ±2 . 



(78) 



We conclude observing that the equations for the velocity-independent source function derived in Chapter 14 of LL04 under 
the approximation of complete redistribution on velocity can be recovered, as a particular case, from the equations presented in this 
work. Indeed, they can be obtained starting from Eq. (1331 (or one of its following reformulations for the case of a plane-parallel 
atmosphere), assuming that the source function Sg and the depolarizing rate <5„ do not depend on the velocity, integrating over the 
velocity components v' a and v' b , and averaging over the Maxwellian distribution of the velocity component u c . 



7. The R\ redistribution phase-matrix 

We now analyze the basic equations that underly this physical problem, following an alternative approach. Instead of eliminating the 
"radiation field variables" in order to obtain a set of coupled equations for the velocity-dependent density matrix elements, we can 
eliminate the "density matrix variables" in order to obtain equations which directly involve the Stokes parameters of the radiation 
field. This will bring us to the definition of a suitable redistribution phase-matrix, referred to in the literature as Ri in the case of a 
Maxwellian distribution of velocities. 

Neglecting the magnetic field contribution, Eq. ( fT/Tb can be easily solved for Sg{v,x). One gets 



X (-1) 2 [^ e (vo)]c.f. + Weo efir(vo) 



(79) 



We can now substitute this expression for the velocity-dependent source function into Eq. ( l23l l. giving the emission coefficient for 
the Stokes parameters at frequency v into the direction il. Taking also into account Eqs. (fTTT > and (ETT i. the expression of the emission 
coefficient can be cast in the form 



Si (v, £1) = J y^— p(v, £1) B T (v ) + y-j- j d 3 vf(v ) 6 {v + v ^ - vj j> ^- g P^fiSl, fi', v) Sj |v& + 



where Pf)(£l, Q', v) is the scattering phase matrix (corrected for depolarizing collisions) which is given by 



" n ' O' 



(80) 



KQ 



1 + 



l + e 



W K (J e , /„)(-! ) G 7?(i, «) T K (j, £2') 



(81) 
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with Wk(J(,Ju) = ( w< fji) ■ ^ e now assume that the velocity distribution is Maxwellian, and we perform the integral in d 3 v, 
neglecting the velocity dependence of the depolarizing collision^. Under this hypothesis, the scattering phase matrix does not 

(c) 

depend any longer on v and can thus be simply written as P\. (SI, £2 ). It is convenient to perform an inversion of the two integrals 

(first we perform the integral in d 3 i> and then the one in AO.') and to introduce a right-handed triplet of unit vectors, (111,112,113), 
defined in the following way: Uj is the unit vector along the direction that bisects the angle (0 < < n) formed by £2 and £2'; 
112 is perpendicular to Ui, lying, as Ui, in the plane defined by £2 and £2', and being directed in such a way that it has a positive 
component along £2 . Finally, 113 is defined accordingly. The unit vectors Ui and 112 are given by 

£2 + £2' £2- £2' 

Ul ~ 2 cos(0/2) ' U2 ~ 2 sin(©/2) ' ( ) 

with the inverse formulae 

£2 = cos(0/2) ui + sin(©/2) u 2 , £2' = cos(0/2) m - sin(0/2) u 2 . (83) 
Writing v — v\ Ui + V2 U2 + 1)3 U3, the delta function appearing in Eq. (f80b is satisfied when 

— ( v - y ) = cos(0/2) v\ + sin(0/2) v 2 . (84) 

Taking into account this relation between v\ and v%, the Stokes parameter Sj(vo + vqv ■ SI' /c, £2'), also appearing in Eq. (l80l l. can be 
written as Sj{v', SI'), where the frequency V only depends on v%, being given by 

V = v-2— sin(0/2)u 2 • (85) 
c 

Integrating first in d«3 and then in dfi (this latter integral being performed taking into account the Dirac-delta), one is left with the 
expression 

£ Kv,n) = ^(x)|^ 7 , M (v-vo)Br(vo) + T ^(£^y^ , (a£2')-^ f — e^ 4 **? Sj(V, £2') , (86) 

1+e 1 + f J ™ ^ ,J nAvo cos(0/2) J v t 



where v\ follows from Eq. ( I84l l. being given by 



1 / 



cos(0/2) \vo 



This equation can be written in an alternative form by a change of variable in the second integral, passing from the variable 02 
to the variable V . This can be done by taking into account that (see Eq. (l85l >) 

^ = l - dv', (88) 

Vt 2 Av D sin(0/2) 

and that, starting from Eqs. (l85l l and (l87l i. the quantity (v\ + v 2 ^)jv\ can be written, after some algebra, in the form 

= , 1 , f(v - v ) 2 + (v' - vo) 2 - 2 (v - v () )(v' - Vb) cos0l . (89) 
v\ Avl sin 2 © L J 

Moreover, by introducing the reduced variables £ and f, defined by 

f = ^p. (90) 
Av D Av d 

taking into account that e,(^, £2) = e,(v, £2) dv/d^ = e,(v, £2) Avd, and recalling the definition of the absorption profile <p(g) (see 
Eq. (l52l), Eq. (l86*T l can be written in the following form 

e ^n) = ^(x)|^^)BHvo) + T ^^^|]^ , (ano ;r ^ J de^^- 2 #^'rt e s0,aM (9i) 

or, using the lexicon of 'redistribution functions', 

Bid, si) = kt(x) j tp{0 B T (v ) + ^ / Z [m > n ' s ^'> n ' } J ' (92) 



3 This is a non-trivial approximation. Indeed, if the colliding hydrogen atoms have an average quadratic velocity w, an atom that is moving 
with velocity v "sees", in his rest frame, that the hydrogen atoms have an average quadratic velocity V^ 2 + w 2 . This can lead to an important 
dependence of the quantities S {K) on v. 
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where 

[Ri(& £2, f , a%j = _i^pg>(n, 12') e- (f2+f ' 2 - 2 «' cos0 )/ sin20 . (93) 

This redistribution matrix was first proposed heuristically by Dum ont et al.l (Il977h . We also observe that for a — > 1 transition (for 

(id 

which Wq - W2 - 1), and in the absence of dep olarizing c ollisions (6„ = 0), the redistribution matrix element [/?/]oo corresponds 
to the redistribution function Ri derived by Hummer ( 1962) for the unpolarized case (see his Eqs. (2.21.2) and (2.21.4))Q 



8. Conclusions 

In this paper we have derived the equations for the Non-LTE problem of the 2 nd kind, taking velocity density matrix correlations 
into account. We considered the basic case of a two-level atom with infinitely sharp upper and lower levels, and we derived the 
statistical equilibrium equations for the velocity-space density matrix, neglecting the generalized Boltzmann term (which is a good 
approximation in the outer layers of a stellar atmosphere). Taking the Doppler effect into account, we derived a set of coupled 
equations for the velocity-dependent multipole components of the source function. Such equations show a coupling between the 
atoms at point x, moving with velocity v, and the atoms at point x', moving with velocity v', such that the difference v - v' (the 
relative velocity) is perpendicular to the direction x - x'. This is a clear consequence of the Doppler effect and of the fact that we 
considered a resonant transition between two infinitely sharp atomic levels, an academic case that nevertheless allows to investigate 
in detail the physics of pure Doppler redistribution in scattering polarization. 

Starting from the statistical equilibrium equations for the velocity-space density matrix, and from the radiative transfer equations 
for polarized radiation, w e derived the r edistribution matrix corresponding to this physical problem. This redistribution matrix 
(heuristically proposed by iDumont et al.l ll 9771) provides an alternative, equivalent description of the problem, and re presents the 
generalization to the polarized case of the (angle-dependent) R; redistribution function derived by iHummerl (1 1962b . If, on the 
one hand, the redistribution matrix formalism allows a very simple and intuitive description of redistribution phenomena, on the 
other hand, the velocity-space density matrix formalism provides a very transparent picture of the physics of the atom-photon 
interaction. The above-mentioned correlations between atoms located at different points of the plasma remain actually "hidden" in 
the redistribution matrix formalism. It should also be observed that only the average effect of such correlations is taken into account 
if the approximate "angle-averaged" Rf A redistribution matrix is considered. The velocity-space density matrix formalism has also 
the advantage of being suitable for taking lower-level polarization into account, and for describing multilevel atomic systems. An 
important point to remark is the appearance, in the statistical equilibrium equations for the velocity-space density matrix, of the 
generalized Boltzmann term. This term allows to include the effect of velocity-changing collisions into the problem, and points out 
the limits of applicability of any theoretical approach in which it is neglected. Velocity-changing collisions are a complex and not 
yet deeply investigated physical aspect that, however, may play an important role in the lower layers of stellar atmospheres. We 
point out that the equations that have been derived in Sect. 5 (in particular Eq. (l35t ) are very general, and can be applied to arbitrary 
velocity distributions and to plasma structures of any geometry. 
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